Measurement scheme to detect α relaxation time of glass-forming liquid
Zhao Xing-Yu1, 2, Wang Li-Na1, 2, Yin Hong-Mei1, 2, Zhou Heng-Wei2, Huang Yi-Neng1, 2, †
National Laboratory of Solid State Microstructures, School of Physics, Nanjing University, Nanjing 210093, China
Xinjiang Laboratory of Phase Transitions and Microstructures in Condensed Matters, College of Physical Science and Technology, Yili Normal University, Yining 835000, China

 

† Corresponding author. E-mail: ynhuang@nju.edu.cn

Project supported by the National Natural Science Foundation of China (Grant No. 11664042).

Abstract

A measurement scheme for detecting the α relaxation time (τ) of glass-forming liquid is proposed, which is based on the measured ionic conductivity of the liquid doped with probing ions by low- and middle-frequency dielectric spectroscopy and according to the Nernst–Einstein, Stokes–Einstein, and Maxwell equations. The obtained τ values of glycerol and propylene carbonate by the scheme are consistent with those obtained by traditional dielectric spectroscopy, which confirms its reliability and accuracy. Moreover, the τ of 1,2-propanediol in a larger temperature range is compared with existing data.

1. Introduction

The glass transition is the crossover between the α relaxation time (τ) of the glass-former and the observation time.[1,2] Therefore, the measurement of τ as well as its mechanism is one of the key issues of the glass transition.[310] So far, there has been no widely accepted theory of glass transition,[36,11] which has led to no exact relationship between τ and temperature (T). Based on available experimental data on τ, a few empirical laws have been put forward, such as the Vogel–Fulcher–Tammann law[1215] (VFT) ( , critical-like behavior[16] ( ), the Waterton–Mauro equation[8] (WM) ( , and the Avramov–Milchev equation[17] (AM) ( , where , A, T0, TC, ϕ, B, C, L, and m are all temperature-independent fitting parameters. The τ values of the VFT and critical-like laws diverge at T0 and TC, respectively, which indicates that there is a phase transition at the temperatures,[3,6] but the τ values of WM and AM only tend to be infinite at 0 K, giving no phase transition above 0 K.

Up to now, the non-zero temperature divergence of τ is still under controversy. For example, Hecksher et al.[5] found that there was not enough evidence of divergence by analyzing 42 glass-formers. However, Drozd-Rzoska et al.[4] showed that the τ value of some glass-forming liquids has a critical behavior. Martinez-Garcia et al.[3] also showed the existence of divergence and a relationship between the glass transition and critical phenomena. Therefore, more systematic experimental data are needed to further understand τ. At present, the wide range of τ is mainly measured by the dielectric spectroscopy,[1820] viscosity,[2124] and diffusion coefficient[2225] methods, and each kind of method contains a combination of a set of technologies. The different techniques have different measurement errors; for instance, high-frequency (greater than about 10 MHz) dielectric spectroscopy has larger errors than low- and middle-frequency dielectric spectroscopy. Moreover, the measurement data of the different techniques have different deviations. So, it is necessary and valuable to explore new measurement methods or schemes for detecting the value of τ.

Consider that the relation between the shear relaxation time (nearly the same as τ in small-molecular liquids) and the viscosity (η) is subject to the Maxwell equation,[1,24] where is the high-frequency shear modulus of the liquid; the relation between η and the diffusion coefficient (D) of ions in liquids satisfies the Stokes–Einstein equation,[2628] where a is the effective radius of the ions, and kB the Boltzmann constant; and the relation between D and the conductivity (σ) of ions in liquids obeys the Nernst–Einstein equation,[2830] where n is the number density of ions, and q the electric charge quantity of the ions.

From Eqs. (1)–(3), it is obtained that where , which is a nearly temperature-independent constant.

So in theory, τ can be obtained by calculating the value of σ. In fact, the experimental results of several materials show that the temperature dependence of τ is the same as that of 1/σ for impurity ions in a large range.[3134] However, both the type and the quantity of impurity ions are unknown; hence it is impossible to control the amplitude of σ and the process of σ will be affected by the α relaxation when σ is small.

In this paper, a measurement scheme is designed to detect the τ of glass-forming liquids by doping an appropriate quantity of given probing ions into the liquids and by using Eq. (4). The τ value of glycerol and propylene carbonate obtained by the scheme is consistent with that obtained by traditional methods. Moreover, the τ value of 1,2-propanediol in a larger time range is given.

2. Measurement scheme

The scheme for detecting the τ value of glass-forming liquids covers three steps as follows.

(i) Dope probing ions (ionic materials, e.g., electrolytes) into glass-forming liquids, and the ion concentration should be low enough, which has no obvious effect on the α relaxation of the liquid.

(ii) Measure the ionic conductivity (σ) of the liquid by low- and middle-frequency dielectric spectroscopy (LMF-DS), which has a high measurement accuracy. Specifically, the imaginary part ( ) of the complex permittivity of the liquid with ions can be expressed as[3537] In other words, where is the imaginary part of the dipole orientation polarization, is the permittivity of vacuum, and ω is the angular frequency of a sine alternating external electric field. As shown in Section 3 (Fig. 5 and 6), in the region where the conductivity plays a dominant role[37] and the electrode effect could be neglected.[38]

(iii) Calculate τ from σ according to Eq. (4), in which c is a reduced constant and is determined on condition that the translation and rotation relaxation times are the same at a certain temperature.

3. Experiments

The typical small-molecule glass-forming liquids, the ionic probe materials, and the concentration (x) values of the ionic probe in the liquids are shown in Table 1.

Table 1.

Glass-forming liquids, ionic probe materials, and concentration (x) values.

.

The parallel-plate capacitor method is used to measure by a Novocontrol Beta-NB impedance analyzer. The measurement frequency (f) values are 1, 10, 100, 103, 104, 105, and 106 Hz. The cooling rate is 1.5 K/min. The values of samples at high frequencies (1 MHz–1 GHz) near the room temperature are obtained by a Keysight E4991B impedance analyzer.

In Table 1, x is the mass percent. The original samples are produced by Aladdin Industrial Corporation (a) and Tianjin Fuchen Chemical Reagent Factory (b).

4. Results and discussion

Figure 14 show the of GLY, PC, and PDO doped with probing ions versus T at a series of f values on cooling. It can be seen that there are two processes in each sample, marked Pα and Ld. The Pα peaks reflect the α relaxation of the liquid, which corresponds to the molecular orientation movement, while the process marked Ld corresponds to the ion diffusion in the liquid. in Figs. 14 is the capacitance of the parallel-plate capacitor without liquid.

Fig. 1. Plots of of GLY doped with 0.1 and 0.3 KOH versus T on cooling.
Fig. 2. Plots of of GLY doped with 1 and 3 NaCl versus T on cooling.
Fig. 3. Plots of of PC doped with 0.2 and 0.6 NaI versus T on cooling.
Fig. 4. Plots of of PDO doped with 0.1 and 0.3 MgCl2 versus T on cooling.

The insets of Figs. 14 show plots of versus T of GLY, PC, and PDO with different x values at f = 103 Hz, and the data at the other f values have the same characteristics. Moreover, all data of each liquid with different x values at the same f coincide well with each other at temperatures both below and near the temperature (Tp) of the maximum of the α relaxation peak, which indicates that the probing ions doped have almost no effects on α relaxation, but obvious differences at temperatures well above Tp. The value of is larger when x is larger, which shows that the Ld process is caused by the movement of the probing ions in the liquid.

Based on the data in Figs. 14, plots of versus T are obtained (Figs. 5 and 6). According to Eq. (6), it can be seen that the overlapping part of in the figures originates from the ionic conductivity,[37] and in the higher-temperature range, the deviations at lower frequencies (curve 1 in Fig. 5(a), curves 1–2 in Figs. 5(b) and 5(c), and curves 1–3 in Fig. 5(d), as well as curves 1–3 in Fig. 6(a), curves 1–2 in Fig. 6(b), and curve 1 in Figs. 6(c) and 6(d)) are due to the electrode effect,[38] i.e., the charge accumulation near electrodes in a period when the probing ions move rapidly at high temperature. Whether the electrode effect is obvious depends on the speed of ion motion and the period ( ) of the external electric field.[38] Only when the period is much longer than the time required for a large number of ions to move to the vicinity of the plate does not satisfy Eq. (7) distinctly. At high temperatures, the electrode effect of low-frequency may be obvious, but when the temperature is lower, the movement of ions will become slower and the electrode effect will be very weak. So low-frequency can also be used to calculate σ as long as it can coincide with those of other frequencies at a certain temperature, i.e., satisfies Eq. (7). From the overlapping parts of , the σ data of the liquids, shown in the insets of Figs. 5 and 6, are given according to Eq. (6).

Fig. 5. Plots of versus T of GLY doped with ions, with inset showing plot of σ versus T.
Fig. 6. Plots of versus T of PC and PDO doped with ions, with inset showing plot of σ versus T.

The solid symbols in Fig. 7 are the relaxation time ( ) of molecule rotations in the liquid (calculated from , and the lines in Fig. 7 represent the τ of GLY, PC, and PDO in a series of x, with c in Eq. (4) chosen under the condition of at low temperature. The data of GLY (295 K) and PDO (299 K) (solid squares in Fig. 7) are obtained by the radio-frequency impedance analyzer. For GLY and PC, the data on from Ref. [18] are also given in Fig. 7 as empty symbols. From Figs. 7(a) and 7(b), it can be seen that the τ value obtained from the ionic conductivity is consistent with the and τ values in Ref. [18], independently of the x value and the types of probing ions, which shows that the scheme proposed here is reliable and accurate. By this scheme, the τ value of PDO in a larger time range than those in Refs. [3944] is given and it is coincident with the value near room temperature, thus further verifying the feasibility of this scheme.

Fig. 7. Plots of τ and versus T of (a) GLY, (b) PC, and (c) PDO.

It is worth pointing out that the measurement scheme proposed in this paper can detect τ values in a larger time range (10−11 s–10−4 s) only by the LMF-DS technique and ionic conductivity mechanism. By contrast, to measure τ values in the same range, dielectric spectroscopy requires the combination of frequency response analysis (∼104 s–10−7 s) (or an AC bridge, ∼10−2 s–10−8 s), coaxial transmission and reflection (∼10−7 s–10−10 s), and quasi-optics (∼10−10 s–10−12 s).[18] The viscosity method needs a combination of rotating (including a few apparatuses, ∼10−13 s–10−5 s), capillary-flow (∼10−13 s–10−8 s), and parallel-plate (∼10−7 s–10−1 s) viscometers.[22,4547] And the diffusion coefficient method requires the combination of diaphragm cell, capillary cell, and refractive index techniques. In these combination measurements, there are some deviations between the data caused by the different techniques or apparatuses. Therefore, the measurement scheme proposed here will be a useful supplement to the above methods.

In the theoretical part of this scheme, equation (1) is derived in Newtonian fluid,[1] and the experimental results are also well consistent with the theoretical results for small-molecule glass-forming liquids.[1,24] Equation (2) requires the diffused particles to be nearly spherical,[24] and equation (3) needs the conductivity to originate from ion diffusion.[48,49] Therefore, in the scheme, both the molecules and the doped ions in the measured liquid will be spherical as far as possible. Moreover, we should point out that a and will significantly change if a liquid–liquid phase transition occurs in the system, i.e., the interaction between molecules varies obviously. In this case, whether the scheme in this paper is applicable still needs further investigating. Of course, the scheme in this paper also needs more experimental work to further ascertain applicability.

From Eq. (7), the measurement error of this scheme mainly depends on the electrode effect and α relaxation. Compared with in Eq. (5), the weaker the electrode effect and α relaxation are, the smaller the error will be. In the experimental analysis, the error is mainly reflected in the coincidence degree of of each frequency used to calculate the σ value. The error of the final experimental data is also related to the measurement error of the instrument.

5. Conclusions

According to the Nernst–Einstein, Stokes–Einstein, and Maxwell equations, a measurement scheme for detecting the α relaxation time (τ) of glass-forming liquids is proposed, which is based on the measured ionic conductivity of the liquid doped with probing ions by LMF-DS. The obtained τ values of glycerol and propylene carbonate by the scheme are consistent with those obtained by traditional dielectric spectroscopy, which confirms its reliability and accuracy. Moreover, the τ value of 1,2-propanediol in a larger time range than that of existing data is given. The scheme can measure τ values in a larger time range (10−11 s–10−4 s) only based on LMF-DS and the ionic conductivity mechanism. By contrast, to measure τ values in the same range, traditional methods require the combination of a few techniques, and there are some deviations between the data of the different techniques. Therefore, the measurement scheme proposed here will be a useful supplement to existing methods.

Reference
[1] Dyre J C 2006 Rev. Mod. Phys. 78 953
[2] Angell C A 1995 Science 267 1924
[3] Martinez-Garcia J C Rzoska S J Drozd-Rzoska A Martinez-Garcia J 2013 Nat. Commun. 4 1823
[4] Drozd-Rzoska A Rzoska S J Paluch M 2008 J. Chem. Phys. 129 184509
[5] Hecksher T Nielsen A I Olsen N B Dyre J C 2008 Nat. Phys. 4 737
[6] Mckenna G B 2008 Nature Phys. 4 673
[7] Novikov V N 2016 Chem. Phys. Lett. 659 133
[8] Mauro J C Yue Y Ellisona A J Guptac P K Allana D C 2009 Proc. Natl. Acad Sci. USA 106 19780
[9] Tan R R Shen X Hu L Zhang F S 2012 Chin. Phys. B 21 86402
[10] Zhang F Chen Y M Wang R P Shen X Wang J Q Xu T F 2019 Chin. Phys. B 28 47802
[11] Angell C A Ngai K L Mckenna G B Mcmillan P F Martin S W 2000 J. Appl. Phys. 88 3113
[12] Scherer G W 1992 J. Am. Ceram. Soc. 75 1060
[13] Vogel H 1921 Phys. Zeit 22 645
[14] Fulcher G S 1925 J. Am. Ceram. Soc. 8 339
[15] Tammann G Hesse W 1926 Z. Anorg. Allg. Chem. 156 245
[16] Colby R H 2000 Phys. Rev. E 61 1783
[17] Milchev I Avramov A 1988 J. Non-Cryst. Solids 104 253
[18] Lunkenheimer P Schneider U Br R Loidl A 2000 Contemp. Phys. 41 15
[19] Kremer F 2002 J. Non-Cryst Solids 305 1
[20] Lunkenheimer P Kastner S Köhler M Loidl A 2010 Phys. Rev. E 81 51504
[21] Smedskjaer M M Mauro J C Yue Y 2009 J. Chem. Phys. 131 244514
[22] Johnson W L 1999 MRS Bull. 24 42
[23] Bocker C Avramov I Rüssel C 2010 Chem. Phys. 369 96
[24] Masuhr A Waniuk T A Busch R Johnson W L 1999 Phys. Rev. Lett. 82 2290
[25] Sangoro J R Iacob C Naumov S Hunger J Rexhausen H R Valiullin V Strehmel Buchner R Kärger J Kremer F 2011 Soft Matter 7 1678
[26] Ngai K L 2011 Relaxation and Diffusion in Complex Systems New York Springer 3 10.1007/978-1-4419-7649-9
[27] Einstein A 1956 Investigations on the theory of Brownian movement New York Dover 12
[28] Einstein A 1905 Ann. Phys. (Berlin) 17 549
[29] Claisse F Koenig H P 1956 Acta Metall. 4 650
[30] Nernst W 1888 Z. Phys. Chem. 2 613
[31] Stickel F Fischer E W Richert R 1996 J. Chem. Phys. 104 2043
[32] Pawlus S Mierzwa M Paluch M Rzoska S J Rol C M 2010 J. Phys.: Condens. Matter 22 235101
[33] Stickel F Kremer F Fischer E W 1993 Phys. A 201 318
[34] Iacob C Sangoro J R Serghei A Naumov S Korth Y Kärger J Friedrich C Kremer F 2008 J. Chem. Phys. 129 234511
[35] Havriliak S Negami S 1967 Polymer 8 161
[36] Richert R 2010 J. Chem. Phys. 133 74502
[37] Wang J Zhao K S Wu L X 2014 J. Chem. Phys. 141 54502
[38] Ishai P B Talary M S Caduff A Levy E Feldman Y 2013 Meas. Sci. Technol. 24 102001
[39] Nielsen A I Christensen T Jakobsen B Niss K Olsen N B Richert R Dyre J C 2009 J. Chem. Phys. 130 154508
[40] Bergman R Svanberg C Andersson D Brodin A Torell L M 1998 J. Non-Cryst. Solids 235 225
[41] Köhler M Lunkenheimer P Goncharov Y Wehn R Loidl A 2010 J. Non-Cryst. Solids 356 529
[42] León C Ngai K L Rol A C M 1999 J. Chem. Phys. 110 11585
[43] Maggi C Jakobsen B Christensen T Olsen N B Dyre J C 2008 J. Phys. Chem. B 112 16320
[44] Meier R Kahlau R Kruk D Rössler E A 2010 J. Phys. Chem. A 114 7847
[45] Dienes G J Klemm H F 1946 J. Appl. Phys. 17 458
[46] Bakke E Busch R Johnson W L 1996 Mater. Sci. Forum. 225�?27 95
[47] Wilson S J Poole D 1990 Mat. Res. Bull. 25 113
[48] Gao J Zhao Y S Shi S Q Li H 2016 Chin. Phys. B 25 18211
[49] Isard J O 1999 J. Non-Cryst. Solids 246 16